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Abstract. We have examined some basic properties of damped Lya systems(DLAs) by semi- 
analytic model. We assume that DLA hosts are disk galaxies whose mass function is generated 
by Press- Schechter formulism at redshift 3. Star formation and chemical evolution undergo in 
the disc. We select modelled DLAs according to their observational criterion by Monte Carlo 
simulation using random line of sights and disk inclinations. The DLA ages are set to be 1 to 3 
Gyr. By best-fitting the predicted metallicity distribution to the observed ones, we get the effec- 
tive yield for DLAs about 0.25^0 . On the basis of this constrain, we further compared our model 
predictions with observations at redshift 3 in the following items: number density; gas content; 
HI frequency distribution; star formation rate density; relationship between metallicity and HI 
column density. We found that the predicted number density at redshift 3 agree well with the 
observed value, but the gas content SIdla is about 3 times larger than observed since our model 
predicts more DLA systems with higher column density. The frequency distribution at higher 
HI column density is quite consistent with observation while some difference exists at lower 
HI end. The predicted star formation rate density contributed by DLAs is consistent with the 
most recent observations. Meanwhile, the connection between DLAs and Lyman Break galax- 
ies(LBGs) is discussed by comparing their UV luminosity functions which shows that the DLAs 
host galaxies are much fainter than LBGs. However, there is a discrepancy between model pre- 
diction and observation in the correlation between metallicity and HI column density for DLAs. 
Further investigations are needed for the star formation mode at high redshift environments. 
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1. Introduction 

Quasar absorption line systems are one of the best objects to trace the physical na- 
ture of the evolving universe. Among various absorption systems, damped Ly a systems 
(DLAs) are the intercepting objects that have highest HI column density. It is a common 
knowledge that DLAs are the progenitors of present-day galaxies. But substantial debate 
continues over exactly what populations of galaxies are responsible for them. Thanks to 
the large facilities available during the past decade, more and more data about DLAs 
are being provided by various observers (see Prochaska et al. 2003). 

In order to have real understanding of the physical nature of DLAs, theorists have been 
using different ways to explain the observed properties, such as simple galaxy model(Hou, 
Boissier & Prantzos 2001; Calura, Matteucci & Vladilo 2003; Boissier, Peroux & Pettini 
2003: Lanfranchi & Friaca 2003); semi-analytical models(SAMs)(Cen et al. 2003; ; Cora 
et al. 2003; Hou et al. 2005; Okoshi & Nagashima 2005); and numerical simulations 
(Gardner et al. 2001; Nagamine, Springel & Hernquist 2003; Churches et al. 2004). 

In this paper, we will adopt a SAM to examine in detail the observed metallicity, HI 
column density and star formation properties of DLAs in the context of standard hier- 



1 



2 



Hou ct al. 



archical picture of structure formation of the universe. The disk galaxy formation model 
with single disks is adopted because we mainly concentrate on HI column densities and 
the cosmic star formation rate density contributed by DLAs rather than their kinemat- 
ics. As an illustration, DLA properties are assumed at redshift z = 3 and the standard 
ACDM cosmogony is adopted. 

2. Model 

2.1. Galaxy formation 

The galaxy formation model in this paper comes from that for disk galaxies suggested 
by Mo, Mao & White (1998, hereafter MMW). In the model of MMW, the halo mass 
function at any redshift z can be described by the Press-Schechter formalism (Press 
&Schschter 1974). The relation between halo mass M and its circular velocity Vc is 
given by M = Vq/(10GH(z)), where G is the gravitational constant, H(z) is Hubble 
constant at redshift z. Details about the model description can be found in MMW. 

The distribution function of halo spin parameter can be described by a log-normal 
function based on numerical simulation. Disks are assumed to have exponential surface 
profiles S(i?) = So exp(— R/Rj), where So and Rd are the central surface density and the 
scale length. The disk global properties can be uniquely determined by parameters A, Vc, 
rrid and the adopted cosmogony, where A is the halo dimensionless spin parameter, Vc 
is the halo circular velocity, irid is the mass ratio of disk to halo. Here, we have adopted 
the nid as a function of Vc in order to consider the galactic winds and outflow(Shu et al. 

2003) . 

After knowing the distributions of Vc and A for halos, we are able to generate a sample 
of disk galaxies by Monte-Carlo simulation in the Vfe-A plane at redshift z ~ 3, which is 
the parent sample for our follow-up DLA simulation. 

2.2. Star formation, chemical evolution and DLA modelling 

The adopted star formation prescription comes from disk galaxy modelling of Boissier & 
Prantzos (2000). Under the approximation of instantaneous recycling, chemical evolution 
in disks can be expressed by the simple closed-box model with metallicity Z to be Z—Zi = 
—phifi, where Zi is the initial metallicity of gas and is assumed to be O.OIZq, p is the 
effective yield, and /z is the gas fraction. We assume that at the initial time (t = 0, z = 3), 
the gas surface density S g o(-R) — So cxp(— R/Rd). Star formation proceeds within disks 
in a typical time scale 1~ 3Gyr (Lanfranchi & Friaca 2003; Dessauges-Zavadsky et al. 

2004) . The effective yield p is assumed to be constant and is obtained by comparing 
the metallicity distributions between model predictions and observations for DLAs. The 
modelled DLAs are selected over the sampled galaxies by random sightlincs penetrating 
disks according to the observed selection criterion, i.e., Ahi ^ 10 2a3 cm~ 2 . Here, random 
inclinations for disks in the sky are considered. 

3. Model results and comparison with observations 

3.1. Observations 

DLAs have shown many observational properties, including metallicities, column densi- 
ties, kinematics, etc. The observed metallicities of DLAs adopted in the present paper 
mainly come from Hou, Boissier & Prantzos (2001) and Kulkarni & Fall (2002) for the 
Zn element. All the data presented by those authors arc compiled from the results of 
various observers. The observed data of HI column densities come from the survey of 
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Figure 1. Distribution of impact parameter(left) and met allicity( right) of modelled DLA popu- 
lation. In the right panel, the solid, dashed and dotted histograms denoting the model prediction, 
observed DLAs with z > 2 and all DLAs, respectively. 



Storrie-Lombardi & Wolfe (2000) (hereafter SW00). Since our model focuses on DLAs 
at redshift z ~ 3, the observed properties of DLAs with redshift z > 2 are selected for 
comparison. 

3.2. DLA host galaxies and metallicity distribution 

In Fig. 1(a), we show the resulted distribution of the impact parameters of selected DLAs. 
It is found that the peak is around 3kpc, which is resulted from the huge amount of small 
halos in PS formalism and the finite radius for individual galaxies that can produce 
DLAs. Because the DLAs are dominated by small galaxies which are always faint, the 
peak implies that the host galaxies of DLAs are difficult to be observed photometrically. 

The only free parameter in the present model is the effective stellar yield p, which 
is determined by comparing the predicted metallicity distribution with observed one at 
redshift greater than 2. The model prediction for the metallicity distributions of DLAs is 
plotted as a solid histogram in Fig. 1(b) while the observed distributions of DLAs with 
z > 2 and all DLAs are plotted as dashed and dotted histograms, respectively. We get 
the effect yield y = 0.25Z Q for the best-fit result with the assumption of star formation 
timescale for DLAs being random between 1 and 3Gyr. 

3.3. SFR density and luminosity function of DLAs 

Based on the selected DLA sample, we can get the predicted SFR density contributed 
by DLAs at z ~ 3 and show it as a cross in the left panel of Fig. 2 which displays the 
cosmic SFR density as a function of redshift resulted from different observations. The 
SFR density for DLAs with z > 2 based on the CII* absorption lines are also plotted 
in the figure as triangles (Wolfe Gawsier & Prochaska 2003, hereafter WGP03). It can 
be found that the model prediction is consistent with observations and supports the 
"consensus" model described by WGP03. 

In the right panel of Fig. 2, we show, with solid histogram, the predicted UV luminosity 
function of selected DLA hosts galaxies. The observed and dust-corrected UV luminosity 
functions of LBGs are also plotted as full circles (data from Aldelberger & Steidel 2000). 
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Figure 2. Cosmic SFR density as a function z with the model prediction contributed by DLAs 
at z ~ 3 as a cross (left panel). The open circles and filled triangles with error bars are the 
observational results from galaxy and from DLAs (WGP03). The right panel shows the UV 
luminosity functions with the solid histogram denoting the modelled DLAs, the dashed and 
dotted lines denote LBGs with and without dust-correction which are taken from Aldelberger 
& Steidel (2000). 

It can be found that the typical -Rab magnitude of the predicted DLA hosts is ~ 30, 
which are much fainter than LBGs with typical Rab ~ 25. This implies that a typical 
DLA host galaxy in our model has its SFR about 100 times smaller than a typical LBG. 
Because the number density of DLA host galaxies is about 0.26 which is about 100 times 
larger than the observed comoving number density of LBGs, our predictions of cosmic 
SFR density contributed by DLAs is similar to that of LBGs at z ~ 3. 

3.4. HI column density distribution and relation to metallicity 

The frequency distribution of HI column density f(Njn,z) for DLAs , which is defined 
as the number of absorbers per unit Nm and per unit absorption distance X, is very 
important for understanding galaxy formation and evolution in the universe. 

We plot the model prediction of /(iVm, z) at z ~ 3 in the left panel of Fig. 3 as a solid 
line while the observed frequency distributions for DLAs with z = 2.5 ~ 3.5 from SW00 
is plotted as circles with error bars respectively. It can be found from the figure that the 
predicted distribution agrees well with observations at the high Nm end. But it is smaller 
than observations at low A^hi end with the maximum difference as large as 3tr, i.e., the 
predicted distribution is a bit flatter than observed ones. The similar discrepancy appear 
in the more complicated numerical simulation study of Nagamine, Springel & Hernquist 
(2003) as well. This could be due to the limitation of our simple model or perhaps this 
could be a failing of the AC DM power spectrum which deserves further investigations. 

Another unusual property of DLAs is that there exists an anti-correlation between 
[Zn/H] and HI column densities, which is independent of redshift as noticed by Boisse et 
al.(1998), who claimed that this anti-correlation is mainly due to observational selection 
effects due to dust. In order to test this bias, Ellison et al.(2001) have surveyed a sample 
of DLAs toward radio selected quasars and found no significant difference in the HI 
distribution to those optically selected quasars. Moreover, dust obscuration has also been 
argued by Prochaska & Wolfe (2002), who made a detailed analysis of dust extinction 



Star formation in DLAs 



5 




Figure 3. Left panel is f(Nm) vs Nm for DLAs with the solid line denoting the model prediction 
at z ~ 3, and the circles with error bars denoting the observations for z = 2.5 ~ 3.5 which are 
taken from SWOO. In the right panel, we show the predicted correlation between metallicity 
[X/H] and HI column density Nm for DLAs with crosses and open circles denoting modelled 
DLAs and observations for z > 2, respectively. The long dashed line is [Zn/H] + Log(A r Hi) = 21. 



contained in DLAs. They found that inferred extinction values and apparent magnitudes 
imply dust obscuration plays a relatively minor effect in the DLA analysis at least for 
z > 2. 

We examine the predicted correlation between metallicity and Hf column density for 
the selected DLA sample in Fig.3(right panel) with the observed results of DLAs at z > 2 
by open circles. As expected, model predictions show an opposite trend compared with 
observations. If we apply the proposed bias of [Zn/H] +Log(A r H i) > 21 in Fig. 3 (long 
dashed line), to exclude the points above this line, the difference still exists. This means 
that the suggestions of Boisse et al (1998) did not help a lot in alleviating the discrepancy 
in our model. 

Another possibility could be the inadequacy of the adopted Schmidt type star for- 
mation prescription in the model. Even for nearby galaxies, the physical basis of star 
formation is still poorly known. Observers have shown various empirical prescriptions 
for star formation in spirals (Kcnnicutt 1998; Wong & Blitz 2002), and most of galaxies 
can be fitted by a Schmidt type law. This has been widely applied in models of galaxy 
evolution. In fact, gas surface density (E g ) includes both the contributions of HI (Shi) 
and H2 (Sh 2 ) with Shi being dominant but the correlation between these two is different 
from galaxy to galaxy. Recent observations of star formation regions in nearby galaxies 
done by Wong & Blitz (2002) showed a complex relationship between SFR and Shi- For 
their spiral galaxies sample (biased to molecule-rich galaxies), SFR shows virtually no 
correlation with Shi, suggesting a maximum HI column density around 10 21 cm~ 2 . This 
is very instructive to the star formation history for DLAs, where the observed HI column 
density seems have an upper limit. 
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4. Summary 

We show in this contribution that the Semi-Analytical Model based on the galaxy 
formation and evolution framework of hierarchical structure formation scenario is quite 
successful in understanding the DLAs properties. Our simple model could well reproduce 
most of the observed DLA properties, such as metallicity and column density distribu- 
tions, star formation rate density contributed by DLAs and so on. However, we found 
that as long as the Kennicutt star formation prescription is adopted, model always pre- 
dict a positive correlation between HI column density and metallicity, opposite to that 
observed. Different explanations exist, for example, observational bias and dust obscu- 
ration. We suggest that the observed trend of anti-correlation could most probably be 
physical (see also Schaye 2001). A most promising explanation might be the inadequacy 
of Kennicutt star formation law in the high redshift environment. More observations and 
theoretical investigations are needed to clarify this trend in the future. 
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